Multiomics-Based Feature Extraction and Selection for the Prediction of Lung Cancer Survival.

Lung cancer is a global health challenge, hindered by delayed diagnosis and the disease's complex molecular landscape. Accurate patient survival prediction is critical, motivating the exploration of various -omics datasets using machine learning methods. Leveraging multi-omics data, this study seeks to enhance the accuracy of survival prediction by proposing new feature extraction techniques combined with unbiased feature selection. Two lung adenocarcinoma multi-omics datasets, originating from the TCGA and CPTAC-3 projects, were employed for this purpose, emphasizing gene expression, methylation, and mutations as the most relevant data sources that provide features for the survival prediction models. Additionally, gene set aggregation was shown to be the most effective feature extraction method for mutation and copy number variation data. Using the TCGA dataset, we identified 32 molecular features that allowed the construction of a 2-year survival prediction model with an AUC of 0.839. The selected features were additionally tested on an independent CPTAC-3 dataset, achieving an AUC of 0.815 in nested cross-validation, which confirmed the robustness of the identified features.

Inferring Bladder Cancer Evolution from Mucosal Field Effects by Whole-Organ Spatial Mutational, Proteomic, and Metabolomic Mapping.

Multi-platform mutational, proteomic, and metabolomic spatial mapping was used on the whole-organ scale to identify the molecular evolution of bladder cancer from mucosal field effects. We identified complex proteomic and metabolomic dysregulations in microscopically normal areas of bladder mucosa adjacent to dysplasia and carcinoma in situ . The mutational landscape developed in a background of complex defects of protein homeostasis which included dysregulated nucleocytoplasmic transport, splicesome, ribosome biogenesis, and peroxisome. These changes were combined with altered urothelial differentiation which involved lipid metabolism and protein degradations controlled by PPAR. The complex alterations of proteome were accompanied by dysregulation of gluco-lipid energy-related metabolism. The analysis of mutational landscape identified three types of mutations based on their geographic distribution and variant allele frequencies. The most common were low frequency α mutations restricted to individual mucosal samples. The two other groups of mutations were associated with clonal expansion. The first of this group referred to as ß mutations occurred at low frequencies across the mucosa. The second of this group called γ mutations increased in frequency with disease progression. Modeling of the mutations revealed that carcinogenesis may span nearly 30 years and can be divided into dormant and progressive phases. The α mutations developed gradually in the dormant phase. The progressive phase lasted approximately five years and was signified by the advent of ß mutations, but it was driven by γ mutations which developed during the last 2-3 years of disease progression to invasive cancer. Our study indicates that the understanding of complex alterations involving mucosal microenvironment initiating bladder carcinogenesis can be inferred from the multi-platform whole-organ mapping.

Hematopoietic stem and progenitor cells confer cross-protective trained immunity in mouse models.

Recent studies suggest that infection reprograms hematopoietic stem and progenitor cells (HSPCs) to enhance innate immune responses upon secondary infectious challenge, a process called "trained immunity." However, the specificity and cell types responsible for this response remain poorly defined. We established a model of trained immunity in mice in response to Mycobacterium avium infection. scRNA-seq analysis revealed that HSPCs activate interferon gamma-response genes heterogeneously upon primary challenge, while rare cell populations expand. Macrophages derived from trained HSPCs demonstrated enhanced bacterial killing and metabolism, and a single dose of recombinant interferon gamma exposure was sufficient to induce similar training. Mice transplanted with influenza-trained HSPCs displayed enhanced immunity against M. avium challenge and vice versa, demonstrating cross protection against antigenically distinct pathogens. Together, these results indicate that heterogeneous responses to infection by HSPCs can lead to long-term production of bone marrow derived macrophages with enhanced function and confer cross-protection against alternative pathogens.

Correction: Kurpas et al. Genomic Analysis of SARS-CoV-2 Alpha, Beta and Delta Variants of Concern Uncovers Signatures of Neutral and Non-Neutral Evolution. Viruses 2022, 14, 2375.

Missing Funding [...].

Detection and characterization of constitutive replication origins defined by DNA polymerase epsilon.

BACKGROUND: Despite the process of DNA replication being mechanistically highly conserved, the location of origins of replication (ORI) may vary from one tissue to the next, or between rounds of replication in eukaryotes, suggesting flexibility in the choice of locations to initiate replication. Lists of human ORI therefore vary widely in number and location, and there are currently no methods available to compare them. Here, we propose a method of detection of ORI based on somatic mutation patterns generated by the mutator phenotype of damaged DNA polymerase epsilon (POLE). RESULTS: We report the genome-wide localization of constitutive ORI in POLE-mutated human tumors using whole genome sequencing data. Mutations accumulated after many rounds of replication of unsynchronized dividing cell populations in tumors allow to identify constitutive origins, which we show are shared with high fidelity between individuals and tumor types. Using a Smith-Waterman-like dynamic programming approach, we compared replication origin positions obtained from multiple different methods. The comparison allowed us to define a consensus set of replication origins, identified consistently by multiple ORI detection methods. Many DNA features co-localized with the consensus set of ORI, including chromatin loop anchors, G-quadruplexes, S/MARs, and CpGs. Among all features, the H2A.Z histone exhibited the most significant association. CONCLUSIONS: Our results show that mutation-based detection of replication origins is a viable approach to determining their location and associated sequence features.

Genomic Analysis of SARS-CoV-2 Alpha, Beta and Delta Variants of Concern Uncovers Signatures of Neutral and Non-Neutral Evolution.

Due to the emergence of new variants of the SARS-CoV-2 coronavirus, the question of how the viral genomes evolved, leading to the formation of highly infectious strains, becomes particularly important. Three major emergent strains, Alpha, Beta and Delta, characterized by a significant number of missense mutations, provide a natural test field. We accumulated and aligned 4.7 million SARS-CoV-2 genomes from the GISAID database and carried out a comprehensive set of analyses. This collection covers the period until the end of October 2021, i.e., the beginnings of the Omicron variant. First, we explored combinatorial complexity of the genomic variants emerging and their timing, indicating very strong, albeit hidden, selection forces. Our analyses show that the mutations that define variants of concern did not arise gradually but rather co-evolved rapidly, leading to the emergence of the full variant strain. To explore in more detail the evolutionary forces at work, we developed time trajectories of mutations at all 29,903 sites of the SARS-CoV-2 genome, week by week, and stratified them into trends related to (i) point substitutions, (ii) deletions and (iii) non-sequenceable regions. We focused on classifying the genetic forces active at different ranges of the mutational spectrum. We observed the agreement of the lowest-frequency mutation spectrum with the Griffiths-Tavaré theory, under the Infinite Sites Model and neutrality. If we widen the frequency range, we observe the site frequency spectra much more consistently with the Tung-Durrett model assuming clone competition and selection. The coefficients of the fitting model indicate the possibility of selection acting to promote gradual growth slowdown, as observed in the history of the variants of concern. These results add up to a model of genomic evolution, which partly fits into the classical drift barrier ideas. Certain observations, such as mutation "bands" persistent over the epidemic history, suggest contribution of genetic forces different from mutation, drift and selection, including recombination or other genome transformations. In addition, we show that a "toy" mathematical model can qualitatively reproduce how new variants (clones) stem from rare advantageous driver mutations, and then acquire neutral or disadvantageous passenger mutations which gradually reduce their fitness so they can be then outcompeted by new variants due to other driver mutations.

Small RNA-Seq Reveals Similar miRNA Transcriptome in Children and Young Adults with T-ALL and Indicates miR-143-3p as Novel Candidate Tumor Suppressor in This Leukemia.

We aimed to identify miRNAs and pathways specifically deregulated in adolescent and young adult (AYA) T-ALL patients. Small RNA-seq showed no major differences between AYA and pediatric T-ALL, but it revealed downregulation of miR-143-3p in T-ALL patients. Prediction algorithms identified several known and putative oncogenes targeted by this miRNA, including KRAS, FGF1, and FGF9. Pathway analysis indicated signaling pathways related to cell growth and proliferation, including FGFR signaling and PI3K-AKT signaling, with the majority of genes overrepresented in these pathways being predicted targets of hsa-miR-143-3p. By luciferase reporter assays, we validated direct interactions of this miRNA with KRAS, FGF1 and FGF9. In cell proliferation assays, we showed reduction of cell growth upon miR-143-3p overexpression in two T-ALL cell lines. Our study is the first description of the miRNA transcriptome in AYA T-ALL patients and the first report on tumor suppressor potential of miR-143-3p in T-ALL. Downregulation of this miRNA in T-ALL patients might contribute to enhanced growth and viability of leukemic cells. We also discuss the potential role of miR-143-3p in FGFR signaling. Although this requires more extensive validation, it might be an interesting direction, since FGFR inhibition proved promising in preclinical studies in various cancers.

Multiomics to investigate the mechanisms contributing to repression of PTPRC and SOCS2 in pediatric T-ALL: Focus on miR-363-3p and promoter methylation.

T-cell acute lymphoblastic leukemia (T-ALL) is a heterogeneous and aggressive malignancy arising from T-cell precursors. MiRNAs are implicated in negative regulation of gene expression and when aberrantly expressed contribute to various cancer types, including T-ALL. Previously we demonstrated the oncogenic potential of miR-363-3p overexpression in a subgroup of T-ALL patients. Here, using combined proteomic and transcriptomic approaches, we show that miR-363-3p enhances cell growth of T-ALL in vitro via inhibition of PTPRC and SOCS2, which are implicated in repression of the JAK-STAT pathway. We propose that overexpression of miR-363-3p is a novel mechanism potentially contributing to overactivation of JAK-STAT pathway. Additionally, by combining the transcriptomic and methylation data of T-ALL patients, we show that promoter methylation may also contribute to downregulation of SOCS2 expression and thus potentially to JAK-STAT activation. In conclusion, we highlight aberrant miRNA expression and aberrant promoter methylation as mechanisms, alternative to mutations of JAK-STAT-related genes, which might lead to the upregulation of JAK-dependent signaling in T-ALL.

Erratum: The origin of bladder cancer from mucosal field effects.

[This corrects the article DOI: 10.1016/j.isci.2022.104551.].

The origin of bladder cancer from mucosal field effects.

Whole-organ mapping was used to study molecular changes in the evolution of bladder cancer from field effects. We identified more than 100 dysregulated pathways, involving immunity, differentiation, and transformation, as initiators of carcinogenesis. Dysregulation of interleukins signified the involvement of inflammation in the incipient phases of the process. An aberrant methylation/expression of multiple HOX genes signified dysregulation of the differentiation program. We identified three types of mutations based on their geographic distribution. The most common were mutations restricted to individual mucosal samples that targeted uroprogenitor cells. Two types of mutations were associated with clonal expansion and involved large areas of mucosa. The α mutations occurred at low frequencies while the ß mutations increased in frequency with disease progression. Modeling revealed that bladder carcinogenesis spans 10-15 years and can be divided into dormant and progressive phases. The progressive phase lasted 1-2 years and was driven by ß mutations.

CRISPRi for specific inhibition of miRNA clusters and miRNAs with high sequence homology.

miRNAs form a class of noncoding RNAs, involved in post-transcriptional regulation of gene expression, broadly studied for their involvement in physiological and pathological context. Inhibition of mature miRNA transcripts, commonly used in miRNA loss-of-function experiments, may not be specific in case of miRNAs with high sequence homology, e.g. miRNAs from the same seed family. Phenotypic effects of miRNA repression might be biased by the repression of highly similar miRNAs. Another challenge is simultaneous inhibition of multiple miRNAs encoded within policistronic clusters, potentially co-regulating common biological processes. To elucidate roles of miRNA clusters and miRNAs with high sequence homology, it is of key importance to selectively repress only the miRNAs of interest. Targeting miRNAs on genomic level with CRISPR/dCas9-based methods is an attractive alternative to blocking mature miRNAs. Yet, so far no clear guidelines on the design of CRISPR inhibition (CRISPRi) experiments, specifically for miRNA repression, have been proposed. To address this need, here we propose a strategy for effective inhibition of miRNAs and miRNA clusters using CRISPRi. We provide clues on how to approach the challenges in using CRISPR/dCas in miRNA studies, which include prediction of miRNA transcription start sites (TSSs) and the design of single guide RNAs (sgRNAs). The strategy implements three TSS prediction online tools, dedicated specifically for miRNAs: miRStart, FANTOM 5 miRNA atlas, DIANA-miRGen, and CRISPOR tool for sgRNAs design; it includes testing and selection of optimal sgRNAs. We demonstrate that compared to siRNA/shRNA-based miRNA silencing, CRISPRi improves the repression specificity for miRNAs with highly similar sequence and contribute to higher uniformity of the effects of silencing the whole miRNA clusters. This strategy may be adapted for CRISPR-mediated activation (CRISPRa) of miRNA expression.

Involvement of miRNAs in cellular responses to radiation.

PURPOSE: Exposure of living cells to ionizing radiation has different consequences, depending on the dose and cell type. Changes in gene expression at the level of transcription and translation, including those regulated by microRNAs (miRNAs), play a role in the intrinsic radiosensitivity of different cells and define their fate, survival or death. The aim of our work was to examine how ionizing radiation may influence the expression of genes regulated by different miRNAs and miRNA biogenesis. MATERIALS AND METHODS: The work was performed on cultured human melanoma Me45 cells, transiently transfected with plasmids containing Renilla luciferase reporter gene targeted by miRNAs Let-7, miR-21 or miR-24. The levels of reporter mRNAs and mRNAs coding for proteins participating in miRNA biogenesis were assayed at different time points in irradiated and non-irradiated cells using RT-qPCR, and reporter protein by luciferase activity assays. MiRNA-targeted motifs in mRNAs coding for proteins engaged in miRNA biogenesis were extracted from the miRTarBase database. RESULTS: Messenger RNA and protein levels of transfected luciferase genes fluctuated in time in patterns that depended on the type of miRNA regulation and changed upon irradiation of the cells. The average levels of reporter mRNAs were higher in irradiated cells, whereas the levels of proteins changed in either direction. Radiation also influenced the levels of miRNAs and the expression of genes engaged in their biogenesis suggesting that the changes in gene expression following ionizing radiation result mainly from these changes in expression of genes regulating miRNA biogenesis and the influence of miRNA on mRNA translation. CONCLUSIONS: Currently, the responses of cells to ionizing radiation are mainly ascribed to changes in their redox conditions and increased intracellular levels of ROS, but the experiments described here suggest that a further important factor is modulation of translation through changes in biogenesis and levels of miRNAs.

In vitro differentiation of ciliated cells in ALI-cultured human airway epithelium - The framework for functional studies on airway differentiation in ciliopathies.

Primary cultures of the human airway epithelium (AE) cells are an indispensable tool in studies of pathophysiology of genetic and environmental pulmonary diseases, including cystic fibrosis (CF), primary ciliary dyskinesia (PCD) and chronic obstructive pulmonary disease (COPD). Air-liquid interface (ALI) culture is the best method to follow the differentiation of ciliated cells, whose dysfunction forms the basis of PCD. Here, we used custom-designed Taqman Low Density Array (TLDA), qRT-PCR-based assay, to analyze expression of 14 AE genes in cells from healthy donors, cultured in ALI settings using Pneumacult medium, with the focus on genes involved in cilia differentiation and in PCD pathogenesis. The results of TLDA assay were compared with the bulk RNAseq analysis, and placed in the cellular context using immunofluorescent staining (IF) of ALI cultured cells. Expression analysis revealed culture time-related upregulation of the majority of cilia-related genes, followed by the appearance of respective protein signals visualized by IF. Strong correlation of TLDA with RNAseq results indicated that TLDA assay is a reliable and scalable approach to analyze expression of selected genes specific for different AE cell types. Characterization of temporal and inter-donor changes in the expression of these genes, performed in healthy donors and in well-defined ALI/Pnemacult culture conditions, provides a useful reference relevant for a broad spectrum of functional studies where the in vitro AE differentiation is in focus.

RNA-seq library preparation for comprehensive transcriptome analysis in cancer cells: The impact of insert size.

With long reads and high coverage, RNA-seq enables comprehensive transcriptome analysis of cancer cells, provided that optimal length of libraries (and their inserts) is assured, to avoid overlap of paired reads and consequent loss of sequencing data. We assessed TruSeq Stranded library preparation protocols (poly(A) enrichment-PA and rRNA depletion-RD) for the thoroughness of transcriptome analysis of a heterogeneous cancer, acute lymphoblastic leukemia. We applied 2x150PE sequencing, >150 M reads/sample on Illumina NovaSeq6000. We show that PA outperforms RD for the analysis of gene expression and structural aberrations. RD is more suitable for detection of various classes of RNAs, mutations or polymorphisms. We demonstrate that reduced RNA fragmentation time (generating longer inserts) positively affects detection of structural RNA changes, without introducing bias into gene expression analysis. We recommend this modification for all RNA-seq studies utilizing reads longer than 75 nt, aimed to go beyond gene expression analysis and to detect also structural changes.

Heat shock factor 1 (HSF1) cooperates with estrogen receptor α (ERα) in the regulation of estrogen action in breast cancer cells.

Heat shock factor 1 (HSF1), a key regulator of transcriptional responses to proteotoxic stress, was linked to estrogen (E2) signaling through estrogen receptor α (ERα). We found that an HSF1 deficiency may decrease ERα level, attenuate the mitogenic action of E2, counteract E2-stimulated cell scattering, and reduce adhesion to collagens and cell motility in ER-positive breast cancer cells. The stimulatory effect of E2 on the transcriptome is largely weaker in HSF1-deficient cells, in part due to the higher basal expression of E2-dependent genes, which correlates with the enhanced binding of unliganded ERα to chromatin in such cells. HSF1 and ERα can cooperate directly in E2-stimulated regulation of transcription, and HSF1 potentiates the action of ERα through a mechanism involving chromatin reorganization. Furthermore, HSF1 deficiency may increase the sensitivity to hormonal therapy (4-hydroxytamoxifen) or CDK4/6 inhibitors (palbociclib). Analyses of data from The Cancer Genome Atlas database indicate that HSF1 increases the transcriptome disparity in ER-positive breast cancer and can enhance the genomic action of ERα. Moreover, only in ER-positive cancers an elevated HSF1 level is associated with metastatic disease.

About 70% of breast cancers rely on supplies of a hormone called estrogen ­ which is the main hormone responsible for female physical characteristics ­ to grow. Breast cancer cells that are sensitive to estrogen possess proteins known as estrogen receptors and are classified as estrogen-receptor positive. When estrogen interacts with its receptor in a cancer cell, it stimulates the cell to grow and migrate to other parts of the body. Therefore, therapies that decrease the amount of estrogen the body produces, or inhibit the receptor itself, are widely used to treat patients with estrogen receptor-positive breast cancers. When estrogen interacts with an estrogen receptor known as ERα it can also activate a protein called HSF1, which helps cells to survive under stress. In turn, HSF1 regulates several other proteins that are necessary for ERα and other estrogen receptors to work properly. Previous studies have suggested that high levels of HSF1 may worsen the outcomes for patients with estrogen receptor-positive breast cancers, but it remains unclear how HSF1 acts in breast cancer cells. Vydra, Janus, Kus et al. used genetics and bioinformatics approaches to study HSF1 in human breast cancer cells. The experiments revealed that breast cancer cells with lower levels of HSF1 also had lower levels of ERα and responded less well to estrogen than cells with higher levels of HSF1. Further experiments suggested that in the absence of estrogen, HSF1 helps to keep ERα inactive. However, when estrogen is present, HSF1 cooperates with ERα and enhances its activity to help cells grow and migrate. Vydra, Janus, Kus et al. also found that cells with higher levels of HSF1 were less sensitive to two drug therapies that are commonly used to treat estrogen receptor-positive breast cancers. These findings reveal that the effect HSF1 has on ERα activity depends on the presence of estrogen. Therefore, cancer therapies that decrease the amount of estrogen a patient produces may have a different effect on estrogen receptor-positive tumors with high HSF1 levels than tumors with low HSF1 levels.

Predicting Time to Relapse in Acute Myeloid Leukemia through Stochastic Modeling of Minimal Residual Disease Based on Clonality Data.

Event-free and overall survival remain poor for patients with acute myeloid leukemia. Chemoresistant clones contributing to relapse arise from minimal residual disease (MRD) or newly-acquired mutations. However, the dynamics of clones comprising MRD is poorly understood. We developed a predictive stochastic model, based on a multitype age-dependent Markov branching process, to describe how random events in MRD contribute to the heterogeneity in treatment response. We employed training and validation sets of patients who underwent whole genome sequencing and for whom mutant clone frequencies at diagnosis and relapse were available. The disease evolution and treatment outcome are subject to stochastic fluctuations. Estimates of malignant clone growth rates, obtained by model fitting, are consistent with published data. Using the estimates from the training set, we developed a function linking MRD and time of relapse, with MRD inferred from the model fits to clone frequencies and other data. An independent validation set confirmed our model. In a third data set, we fitted the model to data at diagnosis and remission and predicted the time to relapse. As a conclusion, given bone marrow genome at diagnosis and MRD at or past remission, the model can predict time to relapse, and help guide treatment decisions to mitigate relapse.

Chronic infection drives Dnmt3a-loss-of-function clonal hematopoiesis via IFNγ signaling.

Age-related clonal hematopoiesis (CH) is a risk factor for malignancy, cardiovascular disease, and all-cause mortality. Somatic mutations in DNMT3A are drivers of CH, but decades may elapse between the acquisition of a mutation and CH, suggesting that environmental factors contribute to clonal expansion. We tested whether infection provides selective pressure favoring the expansion of Dnmt3a mutant hematopoietic stem cells (HSCs) in mouse chimeras. We created Dnmt3a-mosaic mice by transplanting Dnmt3a-/- and WT HSCs into WT mice and observed the substantial expansion of Dnmt3a-/- HSCs during chronic mycobacterial infection. Injection of recombinant IFNγ alone was sufficient to phenocopy CH by Dnmt3a-/- HSCs upon infection. Transcriptional and epigenetic profiling and functional studies indicate reduced differentiation associated with widespread methylation alterations, and reduced secondary stress-induced apoptosis accounts for Dnmt3a-/- clonal expansion during infection. DNMT3A mutant human HSCs similarly exhibit defective IFNγ-induced differentiation. We thus demonstrate that IFNγ signaling induced during chronic infection can drive DNMT3A-loss-of-function CH.

Interferon Gamma Mediates Hematopoietic Stem Cell Activation and Niche Relocalization through BST2.

During chronic infection, the inflammatory cytokine interferon gamma (IFNγ) damages hematopoietic stem cells (HSCs) by disrupting quiescence and promoting excessive terminal differentiation. However, the mechanism by which IFNγ hinders HSC quiescence remains undefined. Using intravital 3-dimensional microscopy, we find that IFNγ disrupts the normally close interaction between HSCs and CXCL12-abundant reticular (CAR) cells in the HSC niche. IFNγ stimulation increases expression of the cell surface protein BST2, which we find is required for IFNγ-dependent HSC relocalization and activation. IFNγ stimulation of HSCs increases their E-selectin binding by BST2 and homing to the bone marrow, which depends on E-selectin binding. Upon chronic infection, HSCs from mice lacking BST2 are more quiescent and more resistant to depletion than HSCs from wild-type mice. Overall, this study defines a critical mechanism by which IFNγ promotes niche relocalization and activation in response to inflammatory stimulation and identifies BST2 as a key regulator of HSC quiescence. VIDEO ABSTRACT.

Discrimination between human populations using a small number of differentially methylated CpG sites: a preliminary study using lymphoblastoid cell lines and peripheral blood samples of European and Chinese origin.

BACKGROUND: Epigenetics is one of the factors shaping natural variability observed among human populations. A small proportion of heritable inter-population differences are observed in the context of both the genome-wide methylation level and the methylation status of individual CpG sites. It has been demonstrated that a limited number of carefully selected differentially methylated sites may allow discrimination between main human populations. However, most of the few published results have been performed exclusively on B-lymphocyte cell lines. RESULTS: The goal of our study was to identify a set of CpG sites sufficient to discriminate between populations of European and Chinese ancestry based on the difference in the DNA methylation profile not only in cell lines but also in primary cell samples. The preliminary selection of CpG sites differentially methylated in these two populations (pop-CpGs) was based on the analysis of two groups of commercially available ethnically-specific B-lymphocyte cell lines, performed using Illumina Infinium Human Methylation 450 BeadChip Array. A subset of 10 pop-CpGs characterized by the best differentiating criteria (|Mdiff| > 1, q < 0.05; lack of the confounding genomic features), and 10 additional CpGs in their immediate vicinity, were further tested using pyrosequencing technology in both B-lymphocyte cell lines and in the primary samples of the peripheral blood representing two analyzed populations. To assess the population-discriminating potential of the selected set of CpGs (further referred to as "composite pop (CEU-CHB)-CpG marker"), three classification methods were applied. The predictive ability of the composite 8-site pop (CEU-CHB)-CpG marker was assessed using 10-fold cross-validation method on two independent sets of samples. CONCLUSIONS: Our results showed that less than 10 pop-CpG sites may distinguish populations of European and Chinese ancestry; importantly, this small composite pop-CpG marker performs well in both lymphoblastoid cell lines and in non-homogenous blood samples regardless of a gender.

Circuits Regulating Superoxide and Nitric Oxide Production and Neutralization in Different Cell Types: Expression of Participating Genes and Changes Induced by Ionizing Radiation.

Superoxide radicals, together with nitric oxide (NO), determine the oxidative status of cells, which use different pathways to control their levels in response to stressing conditions. Using gene expression data available in the Cancer Cell Line Encyclopedia and microarray results, we compared the expression of genes engaged in pathways controlling reactive oxygen species and NO production, neutralization, and changes in response to the exposure of cells to ionizing radiation (IR) in human cancer cell lines originating from different tissues. The expression of NADPH oxidases and NO synthases that participate in superoxide radical and NO production was low in all cell types. Superoxide dismutase, glutathione peroxidase, thioredoxin, and peroxiredoxins participating in radical neutralization showed high expression in nearly all cell types. Some enzymes that may indirectly influence superoxide radical and NO levels showed tissue-specific expression and differences in response to IR. Using fluorescence microscopy and specific dyes, we followed the levels and the distribution of superoxide and NO radicals in living melanoma cells at different times after exposure to IR. Directly after irradiation, we observed an increase of superoxide radicals and NO coexistent in the same subcellular locations, suggesting a switch of NO synthase to the production of superoxide radicals.